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1. Introduction 



Domain walls arise in scalar field theories as solutions connecting two isolated vacua 
which are degenerate. Physical examples can range from a system of liquid crystals [0J to 
defects in cosmological models . A simple way to obtain a theory with degenerate vacua 
is to consider a supersymmetric field theory. In this case supersymmetry guarantees the 
positivity of the scalar potential V(Q), which can be written in terms of superpotential 

2 

W($), i.e. V($) ~ . The location of the minima of the potential are at the 

critical points $ = of the superpotential, such that W'(&k) = 0. 

Starting from the simplest model of a single scalar field theory with a potential in 
1 + 1 dimensions, which allows a single type of domain wall between each of the critical 
points, things get more complicated when we consider theories with multiple scalar fields 
and multiple critical points. In cases where there are more than two degenerate vacua, one 
might consider any pair of vacua and try to connect them with a domain wall (or soliton 
in 1 + 1 dimension). However, this simple-minded construction cannot always be realized 
since there might not always be a BPS solution connecting two given vacua. This can be 
exemplified by the Wess-Zumino(W-Z) model with the following quintic superpotential: 
W($) = $ 5 /5 — <3> 2 /2, which has four critical points, one at $ = and three others 
at vertices of an equilateral triangle. In this theory the domain wall which interpolates 
between $ = and any one of the corners exists, but direct connection of two of the vertices 
does not exist [[|. Therefore such a superpotential only allows for three BPS states and 
not six as one might have expected. (For this particular example, one can actually see 
from surface plot of the potential V(Q>) that there is no BPS path between the vertices of 
the triangle). 

In 1 + 1 dimensions these interpolating BPS solutions are just kinks or solitons. Inte- 
grability conditions for different soliton solutions in 1 + 1 dimensions, interpolating different 
pairs of critical points were studied in Ref. 0, where a soliton which saturates the Bogo- 
mol'nyi bound can best be described as a straight line connecting the critical points in the 
superpotential space, i.e. the W— plane. In fact a very extensive classification program of 
integrable models was carried out in 1 + 1 dimensional theories with M = 2 supersymme- 
try in Ref.[|J. Some of the results there can be used in higher-dimensional theories with 
domain walls because domain walls essentially have one space dimensional dependence, 
which is along the direction separating two domains. One new feature that appears when 
we have more than one spatial dimension is that we can now have intersections or junctions 
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of domain walls 0. We can ask a similar question for the existence of a BPS state between 
critical points each time we encounter a superpotential, and perform an analysis as was 
done extensively in Ref. || . However, it would be desirable to have a more global view in 
the parameter space (i.e. the space of deformations of the superpotential) so that we can 
easily follow the behavior of certain BPS states which are created or destroyed as we move 
around in this parameter space. 

In this paper, we will consider domain walls and their junctions in Af = 1 supersym- 
metric field theories in four dimensions and we analyze under which circumstances certain 
classes of junctions can appear or not. For an appropriate choice of superpotential, such 
domain walls have been shown to arise in the W-Z model and also in SU(N) SUSY QCD 
for which the W-Z model is an effective low-energy theory. Furthermore, it has been shown 
that the W-Z model (at least for a Z3 symmetric configuration of three critical points) 
admits solutions preserving only 1/4 of supersymmetry 0|§||, which were interpreted 
as junctions of three domain walls. More general BPS and non-BPS junctions of the W-Z 
model with a symmetric configuration of critical points where discussed in ||10[| . Re- 



cently nonperturbative junctions of domain wall solutions were also extensively studied in 
SUSY QCD and in the brane world scenarios [I2|[I5|[fH|, where gravitating domain 



wall junctions were considered. 

Another important motivation to study this subject comes from the recent discussions 
of the vacuum and soliton structure of supersymmetric theories in the context of string 
theory compactifications [ 375 j [ |Tfi|] . Consider compactification of Type II, M-, or F-theory on 



some singular noncompact Calabi-Yau n manifold with some background flux of Ramond- 
Ramond field, say G. (For F-theory, we need elliptically fibered Calabi-Yau manifold, and 
in addition we need both NS and RR fluxes.) Nonvanishing R-R flux is needed to cancel 
the tadpole anomaly |17j , while taking a singular limit of a Calabi-Yau manifold leads to a 



decoupling of gravity in the effective field theory in the lower dimension [[Tj| . Domain walls 
are identified with D-branes (or M-branes for M-theory) wrapped on supersymmetric cycles 
and in crossing such a brane the flux (of the appropriate field) jumps, so the different values 
of the flux correspond to different vacua. For supersymmetric vacua certain conditions has 
to be imposed on G |H| . These constraints can be realized by interpreting G as giving rise 



to an effective superpotential of the lower-dimensional theory which is of the form 

W = j A AG, (1.1) 



where A is either the holomorphic n-form O or some appropriate power of the Kahler 
potential /C. i For compactification of Type II, M-theory or F-theory on singular Calabi- 
Yau manifolds this analysis leads in certain cases to an identification of the corresponding 
low-dimensional theories as specific non-trivial conformal field theories, depending on the 



singularity in question. As an example, it was shown [fL6] that Type IIA compactified on 
a Calabi-Yau four-fold with A n singularity gives an Af = 2 Kazama-Suzuki model [^T] 
two dimensions. 



m 



In this paper, we will concentrate on W-Z models in four dimensions (with four su- 
percharges), though much of the analysis can be applied in three and two dimensions as 
well. We analyze the appearance of BPS domain walls and junctions for massive defor- 
mations away from the conformal point. In section 2, we review the possibility of central 
charges of the M = 1 superalgebra in four dimensions and their interpretation in terms of 
domain wall and junction charges and also the BPS condition for the domain walls and 
their junctions. In section 3 we review the derivation of W-Z models in D = 2, 3 from type 
IIA or M-theory and discuss some relations between the geometry of the Calabi-Yau man- 
ifold and the solutions of the BPS equation in lower dimensions. We also comment about 
generating superpotentials in F-theory. In section 4, we collect the rules for the counting 
of BPS states, which are used in section 5 in studying massive deformations of the W-Z 
model with a general quintic superpotential. Finally, section 6 contains our discussions. 



2. Supersymmetry Algebra and the BPS Condition 

We start by recalling the structure of the M = 1 supersymmetry algebra in 3 + 1 
dimensions and how the possibility of domain walls and junctions of domains walls can 
be analyzed directly from this algebra. (For further discussions of the M = 1 algebra in 
£> = 4, see 0@[0|23). 

The M = 1 supersymmetry algebra in D = 4 allows central charges which correspond 
to tensions of BPS domain walls and junctions of them [24|[p2|: 



{Q ai Q l 3} = 2i(a k a°)Je lP Z kl 
{Q a ,Q d } = 2K Q P^ + a^Y fc ), 



2 Note that this is related to the theory of calibrations: A is the calibration and for A = Q, 
these potentials are related to Lagrangian submanifolds and give rise to chiral superfields, while 
if A = K7 they are related to holomorphic curves and lead to "twisted" chiral superfields. 
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where k = 1,2,3 and \x = 0, . . . , 3. The Zk (which are complex charges) have an inter- 
pretation as domain wall charges and Yk (which are real charges) as the junction energy, 
which can be either positive or negative [p2j . 

The relations between the superpotential and the central charges are given by 



where <f> is the scalar component of the chiral superfield, and the Kahler metric is derived 
from the Kahler potential K via = d 2 K j 'dcjfdcffi . The central charges Zk depend 
only on the difference between the values of the superpotential at spatial infinity. If we 
have a single domain wall - which is only nontrivial in one dimension - then Yk vanishes 
for all k and Zj is nonvanishing (for some j) since the Zj central charge depends on the 
spatial derivative in the Xj direction. For junctions of domain walls to be possible, one 
first of all need more than a single (real) scalar field, since else Yk will vanish identically. 
Furthermore, as one can clearly see, Yk is nonvanishing only when the field configuration 
at infinity is nontrivial in two dimensions. If we have two spatial dependences as for a 
domain wall junction solution, then we will generically have two of the Z^'s nonvanishing. 
When the Kahler metric is trivial, Yk is just a surface integral at the infinity. 

When we start from a M = 1 theory in D = 4 we originally have four supercharges. 
Domain wall configurations with nonzero Z^s and vanishing Yk has two conserved su- 
percharges, thus are 1/2 BPS states, whereas when there is nonzero Yk there is only one 
combination of the four supercharges which can survive. This leads to a 1/4 BPS state. 

The BPS equation for a single static domain wall of the W-Z model (dimensionally 
reduced to two dimensions) is given by: 



where the prime denotes derivative with respect to <fi = <fi(x), and x is a coordinate and a 
is an arbitrary phase. A domain wall solution of mass M saturates the bound M > \T\, 
where T is the topological charge associated with the wall and has a = argT. The BPS 
equation for a domain wall junction can be derived in higher space dimension as in |7| 
and is completely analogous to Eq.(|2.4j). In particular, if we suppress spatial dependences 
other than two of them, say x and y, then the BPS equation becomes 




(2.2) 




(2.3) 



d x( p = e ia W 




(d x - i 



d y )</> = e ia W. 
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The BPS solution saturates the bound M > \T\ + Q 1 where Q is the junctions charge. 
When there is only on spatial dependence, e.g. d y 4>(z) = 0, this reduces to Eq. ([2.4]). Note 
that this junction is an object in a three-dimensional theory and not a two-dimensional 
theory as the one discussed in H. 



3. Wess-Zumino Models from Calabi-Yau Compactifications 

The W-Z model we will consider will be a field theory in D = 4 dimensions with 
M = 1 supersymmetry. It has a superpotential W($) which is of the form 

m 

W'($)=Cl[(<5>-\ l ), (3.1) 

i=l 

where C is a constant and Aj, i = 1, . . . , m are the locations of the critical points in the 
field space. The Z m symmetric case, for Aj = |Ao|e 27ri,/m , (i = 1, • • • , m), were considered 
in connection to the M = 1 supersymmetric SU(N) YM theory in the large iV limit 
Although it is believed that W-Z models for m > 3 will flow to trivial IR theories for 
D > 3, it can become relevant as a perturbations to certain fixed points. Furthermore one 
can have non-trivial brane configurations realizing these higher order potentials |f25[ . We 
will comment on this in the last section. 



As discussed in |T^] Wess-Zumino models can arise in Calabi-Yau compactifications 
of M/Type IIA theories. The locations of the isolated singularities correspond to the 
locations of the critical points. In the case of an singularity, the local geometry of the 
Calabi-Yau n-fold is described by an equation of the form 

-P m ( Zl )+z 2 2 + ... z 2 +1 = (3.2) 

with P m (zi) a generic polynomial of degree m = k + 1 in z±: 

Pm(zi)=Y[(z 1 -a i ). (3.3) 

i 

In the above a, are the locations of the singularities. When we fix z\ we can regard |P m (^i)| 
as the radius of the n — 1 sphere which is the nontrivial cycle of the manifold. Since the 
mass of the brane wrapping around the singularity will be proportional to the volume of 
the sphere, and this mass will give also the tension of the domain wall in the lower effective 
field theory, we get the following fact that the superpotential W is related to P m through 



n — 2 



dW = P m 2 dz u (3.4) 
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where the right hand side is basically the volume of the sphere. So, for Calabi-Yau fourfold 
compactifications we recover the W-Z superpotentials. 

Instead of going directly to the case of a quintic superpotential, we would like to 
discuss some general features of the solutions of the BPS equation corresponding to com- 
pactification on a general Calabi-Yau n-fold. BPS states are in this case identified with 
wrapped n-branes in a Calabi-Yau n-fold near an isolated singularity []TB|]. The kind of 
singularities we will be looking at are the Ak singularities, which are described by the 
Eq.( p.2| ). Any such Calabi-Yau manifold has a unique n-form Q which determines the 
volume of a cycle C. The condition for a cycle to be minimal is that its volume saturates 
the inequality 



V 



[ M > 


[ 


Ic 


Jc 



(3.5) 



Now, the problem of minimizing this volume can be mapped to a problem in the complex 
zi-plane as follows. One considers the n-cycle to be an n-sphere S n , which is locally of 
the form S™ -1 x S , i.e. as an 5 <n_1 -sphere fibered over a real curve in the ^i-plane. The 
local volume of this 5 <n_1 -sphere is determined by z\ + • • • %n+i = P m {z\) and so vanishes 
at the roots of P(z±), which are identified with the critical points of the superpotential 
W(z±). With this choice of local coordinates on the singular Calabi-Yau, the expression 
for the holomorphic n-form is 

_ dzi--- dz n 



Z n -\-i 

idz\ ■ ■ ■ dz n (3-6) 



z( + ■ ■ ■ + z 2 n - P m ( Zl ) 

The condition for a cycle to be supersymmetric is that the image of the path is a straight 
line in the flat VF-plane, where W is defined through the relation in Eq. ( p.4j) . This comes 
from minimizing the l.h.s. of ( |3.5| ) with the expression ( |3.6|) for O. The BPS condition is 
then: 

W(z(t))= p—dz = at, (3.7) 

J ZQ 

where t parametrizes the curve connecting the two critical points. To obtain the BPS 
states one should therefore solve the first-order differential equation: 

„ n-2 dZ 

P ^M =a (3 ' 8) 



with the boundary conditions that z(t) should begin and end at the roots of P(z) (or 
rather of dW(z)). Near a root, which we take to be at z = 0, one is solving an equation of 
the form 

for which the solution is 

/ in \ 2/n 

z = {-at) . (3.10) 

In the case of a Calabi-Yau four-fold we see that there are four solutions for any a and 
that the corresponding curves intersect at an angle of 90°. 

Now we will discuss how to construct domain walls in such Calabi-Yau compact- 
ifications and we will follow the discussion in For more details, see also [2B|[E7|. 



Consider compactification of M-theory on some Calabi-Yau four-manifold Y with some 
background flux of the three-form potential C which couples to the membranes (see Table 
1, in which we summarize the construction of vacua and domain walls in compactification 
of M/IIA/.F- theory with background fluxes as in |0J]). These C-field are classified by a 
class £ G H (Y; Z), which in turn defines a lattice T* = H^(Y; Z). This set of data speci- 
fies a choice of vacuum. Now, to make a domain wall in R 3 , one considers a fivebrane with 
worldvolume R 2 x S, where S is a four-cycle. Such four-cycles are classified by H^Y; Z), 
which defines a lattice T dual to F*. So when crossing such a domain wall £ changes - and 
the possible values of £ are classified by F*/F. The effective superpotential is obtained as 
follows. When crossing a domain wall, the change in the superpotential is 



AW = / OA AG, (3.11) 
*x 



2n 



where G = dC. Here X = R 3 x Y. This superpotential will then account for the restric- 
tions on G (which are implied by having vacua with supersymmetry) mentioned in the 
introduction. 

Compactifying Type II A on Y, we have X = R 2 x Y. And to specify a vacuum we 
should also specify the topological class of the (7-field, which is now the RR four-form 
field and takes values in if 4 (Y;Z). To make a domain wall, one now has four-branes 
with worldvolume R x S. Again the possible four-cycles S are classified by H^(Y; Z) and 
therefore £ takes values in T*/r. Compactification of Type IIA on Calabi-Yau four-fold Y, 
with Ak singularity, will then give an effective two-dimensional theory with superpotential 
determined by (fTl) for n = 4. This is precisely the superpotential discussed in the 



following sections, and here we can of course have domain walls between different vacua. 
But we will not have junctions. 

The story for F-theory |2l| compactifications is slightly different. First of all, for F- 
theory compactification on R 4 x Y we need Y to be an elliptically fibered four-manifold. 
The flux $ discussed in [|16]| now has contributions from space-filling threebranes and not 



membranes as in the compactification of M-theory. The analog of the G-field now becomes 
both NS and RR three-form fields, H NS and H R , from the Type IIB theory. F-theory 
on R 4 x Y can be described as Type IIB with certain (p, q i )-sevenbranes on a locus L C B, 
where B is the base of the elliptic fibration. However, in this situation, one can find a 
simpler description: This F-theory compactification with singularity can be reinterpreted 
as Type IIB with a D7-brane with worldvolume R 4 x L, where L is a complex (singular) 
surface inside C 3 (see Table 1). This specifies a choice of vacuum. One has a £7(l)-gauge 
field on the D7-brane and so this vacuum is characterized by the first Chern class, or an 
element f of the lattice T* = H 2 (L; Z). 



M— theory 


IIA 


— theory 


F— theory 


Vacuum : R 3 


x Y 




R 2 : 


< Y 


R 4 


x Y 


G = 


dC 




G = 


dC 


D7 = R 4 


x L 


GeH 4 (Y;Z) = 


J™ 1 * 


G E H 4 (Y 


Z) = 




F 6 H 2 (L;Z) = 




Domain wall : M5 = R 2 


x S 


DA 


= R 


x S 


D5 = R 3 


x V 


[S] E H 4 (Y;Z) 


= r 


[S] 6 H 4 (Y;Z) 


= r 


[8V] E H 2 (L;Z) 


= r 



Table l:M/IIA/F-theory on Calabi-Yau four-fold Y. 



How do we construct domain walls? Take a D5-brane, which can end on the D7-brane, 
with worldvolume R 3 x V, where V is a three-manifold whose boundary should be in 
L (since the D5-brane ends on the D7-brane). This boundary defines a topological class 
[dV] G H2(L; Z), i.e. in the lattice V = H 2 {L; Z) dual to T*. Crossing the domain wall, 
the Chern class changes by the amount [dV]. Again £ takes values in F*/F. We also need 
to specify the local geometry of Y. For elliptic four-fold singularity one has the description 

y 2 = x s + 3ax 2 + H(z 1 ,z 2 ,z 3 ), (3.12) 
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where if is a polynomial in Z\, Z2, Z3. The equation for L then becomes simply H = and 
to describe an A^-singularity one should then choose: 



H = z k+1 + 4 + 4- (3.13) 

It would be desirable to have an explicit computation of the superpotential in F-theory 
generated by the inclusion of iJ-flux and with A — D — E-type singularities. For that 
one could start with Type IIB on Calabi-Yau three-fold as in [29|, where W = J Q A 



(tH ns + H R ), and then lift this construction to F-theory. Note, however, that not all Y 
will generate a nontrivial superpotential 



4. Rules for the Construction 

Now we will discuss the rules for finding the number of BPS states for different values 
of the perturbation parameters, which translates to varying the positions of the critical 
points. 

1) What are we constructing? 

From the BPS equation one can easily show that the BPS solution trajectories are 
straight lines between critical points in the VF-plane ||. However the inverse image of a 
certain straight line - connecting, say W{i) and W(j) - might not lift back to the field 
space as a curve connection the vacua and thus does not correspond to a BPS solution. To 
count the number of actual solutions connecting i and j, one starts with the "wavefront" (or 
sphere) of all possible solutions emanating from i with fixed values of W, denoted by Aj, 
and the same for the critical point j. The number of solutions is then exactly the number 
of points at which Aj and Aj intersect [|J(note that the intersection number depends on 
a choice of orientation and what we really are computing is a weighted sum[^l|). This 



defines the intersection number fiij = Aj o Aj as a quantity which is invariant under 
small perturbations of the superpotential since it is integer. However, as we vary the 
superpotential the critical points will move around in the VF-plane and when a third root 
k crosses the straight line connecting i and j the number of BPS solutions connecting i 
and j can obviously change. Precisely how this number changes can be derived using the 
Picard-Lefschetz theorem and is given by Q 

/4? = Vij =•= Vik ■ Vkj (4.1) 
9 



Here the ± sign depends on the ordering of ikj in the triangle defined by the three roots 
before k was crossing the line between % and j. (Physically, this change in the intersection 
number, as one root crosses the line between two other roots, can be understood in terms 
of the Hanany-Witten effect |3!|] []33[].) In principle on can determine these intersection 



numbers by solving the so-called tt* equation @[|J for a fixed choice of superpotential. 
But in our case we vary the parameters in the superpotential and it is more straightforward 
to look at conditions on masses of BPS solutions (and phases of the topological charges) 
to determine which kind of junctions exist or not. So we are in a certain sense trying to 
give a unified description of the cases considered in Ref. 0. 

So, it would be nice to have the form of functions of the parameters of the 

theory. Since it takes integer values, it is stable under small perturbations and changes 
only by an integer, and the best way to represent it would be to find the boundaries in the 
parameter space where the jump in the values happens. (This is called a separatrix curve.) 
Then we can specify the values of //ij's inside each domain separated by the boundaries 
in the parameter space graphically. In the case of W-Z models we have = or 1. 
Crossing a boundary induces a change in the number of BPS state of ±1. So the graphical 
representation will be as follows. We will denote the critical points as dots. Then we will 
link the critical points i and j by a solid line if = 1. We will not link them if /j,^ 
vanishes. There will be at least one line coming from each critical point. (The connectivity 
is quite analogous to Dynkin diagrams.) So, if there are k critical points, there will be a 
maximum of k(k — l)/2 BPS states and a minimum of k — 1 BPS states since all critical 



points can be connected through a sequence of BPS solutions |34 
2) What determines the separatrix equation? 

Observe that the topological charge associated with two critical points i and j is 

Tij = 2e l ^ {w{z ^- w{zi)) \W{ Zj ) - W{ Zi )\ (4.2) 

and so is a complex number. The mass M of a domain wall is bounded by the absolute 
value of the topological charge T: 

M > |T|, (4.3) 

and is saturated by a solution of the BPS equation. Now consider a situation where i and 
j and also j and k are connected by a domain wall solution with BPS masses My, Mj^ 
and topological charges Ty, Tjk- Let us consider the possibility of a BPS object between 
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i and k. The possible BPS mass of such a solution is always bounded by the following 
simple inequality: 

M lk = \Tij + T jk \ < \T ZJ \ + \T jk \ = Mij + M jk . (4.4) 

The inequality is saturated only when the phases of Tjj and Tj k are the same. When the 
equality ( |4.4j ) is saturated, such that M\ k = Mij+Mjk, then the domain wall with mass Mik 
decays into the two other domain walls. Since the phase of the topological charge comes 
from the argument of the difference of the superpotential, we can calculate the boundaries 
in the deformation parameter space where different solitons are created or destroyed as 
we change the parameters. Each such boundary is determined by three critical points and 
determines whether a solution between a particular pair of them becomes unstable or not. 
The entire parameter space will therefore be divided into many different domains and each 
domain will have the same number of possible BPS solutions. 

3) To map the entire parameter space we pick a point in the space where the BPS 
configuration is easily determined. As we move across a boundary a certain state can be 
created (if it was not there) or destroyed (if it was there). This technique will be applied in 
the next section where we find the separatrix curves for a general quintic superpotential. 



5. Finding the BPS Configurations 

5.1. Quartic Superpotential 

The simplest nontrivial superpotential is of course one with two critical points. This 
allows a single BPS state and hence a single domain wall. Next would be one which has 
three critical points. In this case of k = 3 roots, and actually for all k > 3, one can argue 



that any pair of roots can be connected through a sequence of domain wall solutions [34 



By rescaling and fixing the value of the field we can fix two of the critical points to be, 
say at z\ = — 1 and Z2 = 1. The third critical point can be at an arbitrary point in the 
complex plane, say at Z3 = fi (this case is discussed in detail in [|J). When \x is a real 
number, (j, > 1, the critical points in the W— plane will be colinear and the only straight 
line connecting z\ with 23 will be through z-i- So there can only be two types of domain 
walls. The same conclusion - i.e. that there are only two BPS states - can be drawn when 
\p\ < 1 for real //, and also for (j, < — 1. Let us now see what happens when we move away 
from the real line, holding fixed Z\ and Z2, for the case of — 1 < ji < 1. As ji = n\ +1^2 (pi, 
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H2 are real numbers) moves away from the real line, the number of BPS states stays the 
same until we reach a boundary in the complex \x plane where a new BPS state appears, 
arising from the domain wall between z\ and z 2 ■ This boundary is defined by the condition 
that the phases of the topological charges T i3 and T 32 are the same [Q. When the phases 
are the same then the inequality of the masses saturate and we have M\ 2 = M13 + M32 
(Mij is the mass of the soliton connecting the roots Z{ and Zj). Similar boundaries can be 
found for the initial cases of \Re(fx) \ > 1, determined by the equality of masses of M i2 and 
M 23 or M31 and M 12 . 

There is a reflection symmetry of the boundaries in the real line. These three bound- 
aries together form the separatrix curve and the equation can be written down as the 
following condition on the real and imaginary parts of \i: 

3/4 + 2/^-/4-6/^-6/^ + 3 = 0. (5.1) 

Note that this equation does not distinguish which BPS state melts away as we cross a 
boundary. Different branches of eq. (|5.1| ) will correspond to one of the boundaries which 
we discussed above, obtained from the relations between the possible masses of domain 
walls. So if we put = + Mjk — Afjfcj the separatrix equation will be equivalent 
to i 7 'i23-f r 'i32-f3i2 = 0, after some algebraic manipulations. This observation will be quite 
crucial in identifying various BPS states in the cases with more than three critical points. 
The real line will appear as a solution of the separatrix equation, but it will be a line of 
marginal stability, so the number of BPS states do not change as we cross the real line. 
The connectivity of the roots for the quartic superpotential is therefore very simple: either 
any root is connected to any other root (for a total of three BPS states), or two of the 
roots are not directly connected (for a total of two BPS states). This result is given in 
Fig. 4 of Ref. ||. The connectivity signals a possible BPS state. It can be occupied or be 
vacant. Now when the occupied BPS states are such that we have an enclosed domain, 
then we have a junction of the three domain walls and a 1/4 BPS state. If they do not 
enclose a separate domain, say just two of the edges of a triangle, then we have two BPS 
domain walls which never join and the whole configuration will be 1/2 BPS. So when the 
positions of the critical points are more or less colinear (in W-space) domain wall junctions 
do not develop. This can be used in the cases with more than three critical points, where 
the positions of three particular ones will more or less follow the pattern described above, 
although the very existence of the other critical points do interfere with the detailed shape 
of the separatrix curves. 
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5.2. Quintic Superpotential 

Next we analyze the D = 4 W-Z model with a general quintic superpotential, 



where ctj, i = 1, . . . , 4 are complex deformation parameters. Critical points of the superpo- 
tential are points z a where dW(z a ) = 0. The possible connections of the critical points in 
this case are shown in Figure 1, where each dot represents a root i (or critical point) and 
each line represents a possible domain wall solution interpolating between critical points i 
and j. When such a line exists between two given roots, \fiij\ = 1, and it thus represents a 
possible BPS state. When there is no line between two given roots = and there can 
be no BPS state. Therefore it is easy to see that Figure 1 exhausts all possible connections 
of critical points. So one expects that in some domain the number of BPS states is the 
smallest possible, namely three (as in Figs 1-A,- B), while in some other domain the max- 
imum number of BPS states, namely six (as in Figure 1-F) is obtained, depending on the 
choice of superpotential, i.e. deformation parameters a.i. However, Figure 1-C deserves 
some further comments. In the following we will see that in no finite domain of defor- 
mation parameters will the connectivity be as in C. This is actually easy to understand 
geometrically. In such a four-gon - defined by roots i, j, k and I - are i'j'k' connected for 
any cyclic permutation of the four roots but %' and k' are not connected. So all the angles 
of i'j'k' has to be at least 90°. But the sum of the angles of the 4-gon is 360° and we have 
a contradiction. What about the case of k > 4 number of critical points? For a /c-gon, the 
sum of angles is (k — 2) • 180°. For a configuration with no "internal" BPS solitons the 
sum of angles should be at least k ■ 90°. For k > 5 one might have such domains. 



4 




(5.2) 



i=i 




c) 
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Fig.l: Possible connectivities of four critical points in the case of a quintic superpotential. 

The soliton structure of any such massive deformation of a conformal theory is charac- 
terized by the matrix S = 1 — A, where A is an upper triangular matrix whose elements 
Aij for i < j are exactly A^ = = Aj o Aj ||. However, this matrix does not in itself 
classify the possible junctions. Precisely for this reason will Figure 1-F need some further 
comments, if the actual location of the critical points is as in Figure 1-F (triangle), then 
one can obtain a junction of domain walls by occupying all six states (this junctions will 
look like a circle with three legs coming out). But imagine that the locations of the critical 
points are as in Figure 1-F (square) with the inclusion of the two BPS states connecting 
diagonal corners. Occupying all these states would not give rise to a stable junction. 

In Eq. (|5.2| ) we fix two of the critical points to be at z = ±1, so that the four critical 
points are located at 

zi = -1, z 2 = 1, z 3 = //, z 4 = A, (5.3) 
corresponding to the superpotential which takes the following form: 

W = z 5 - -(/I + X)z 4 + -(/jlX - l)z 3 + -{fi + X)z 2 - hixXz. (5.4) 
4 3 2 

We thus have two complex parameters \x and A to vary, and in general it is not easy to 
visualize different domains in this space of parameters. A systematic way is to fix one of the 
complex parameters, say \± and have a sliced view of the separatrix walls. We will consider 
a few representative values of [i: 1) the case where three points z = —1, z = 1, z = fi are 
at vertices of an equilateral triangle, (This includes the case we already discussed in the 
introduction which corresponds to the situation where the fourth critical point is at the 
center of the triangle. For this case we already know the possible connectivities of the 
critical points and we can use it as the 'initial data' for our analysis.) 2) the case where 
three points are colinear on the real axis and finally 3) the case which includes the Z4 
symmetric case. 

For a generic configuration of roots (i.e. when z% is not colinear with z\ and z%) one 
can obtain the complete set of separatrix curves as follows. Pick any two roots Z{, Zj 
(i > j) and consider the basic separatrix curve joining them as defined by the equation 
F^jF^jiFj^i = 0. Then the condition that the product of all these groups of terms vanishes 
is the equation for the "complete" separatrix curve, just as it is in the case of a single pair 
of roots when we have a quartic superpotential. Now we will focus on the three cases. In 
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the first case we take the three fixed roots to be at the vertices of an equiiateral triangle, 
i.e. Z\ = —1, Z2 = +1, zs = (see Figure 2). 




Fig. 2: Z3-symmetric case. Connectivity of roots depending on the value of Z4 = A. The three 
fixed roots are located at zi = — 1, zi = 1 and Z3 = i\f?>. 

In this case there is a Z3-symmetry generated by rotations of 27r/3 in the center of the 
triangle. In the second case we take the roots to be colinear z\ = — 1,^2 = +1,-23 = +3 
(see Figure 3). In this case there is a Z2 symmetry generated by reflections along the 
vertical line A2 = 0. The last configuration is where z\ = —1, zi = +1 and Z3 = — 1 + 2i 
(see Figure 4) and so contains the Z4-symmetric superpotential (for Z4 = 1 + 2i) which 
has been much studied. 

Before going into details with the different phase diagrams and determining in which 
domains we have how many BPS states and so forth, we start with a global view (i.e. 
far away from the origin). What determines the angles between the curves of marginal 
stability? For that we will take a long-distance view of the separatrix curves. This limit 
corresponds to both Ai and A2 large. For any fixed value Z3 = fx, one can write down the 
separatrix equation as a sixth order equation in Ai and A2 (for example for the pair of 
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roots (z2, Z4) and (24,2:3) as follows): 

= - ^Af A 2 (-l + pf(3 + fi) + |AfA 2 (-l + ^) 3 (1 + 3// + p 2 ) 

- ^A?A 2 (-1 + p) 3 (-5 + p + 3p 2 + p s ) 

- 4a 2 A 2 (-1 + /") 3 (6 + 18// + 6// 2 + X 2 2 (6 + 3// + // 2 )) 

lo 

- ^A 2 (-l + /") 3 (15// 2 (2 + //) + A^)(-3 + 6// + 2// 2 ) + 5A 2 (-2 + 9// + 3// 2 )) 

lo 

+ — AiA 2 (-l + /") 3 (3A2(3 + //) + 3/x(8 + 9// + 3// 2 ) + A 2 (15 + 13// + 9p 2 + 3// 3 )). 

ou 

(5.5) 

The sliced view of this separatrix equation will be shown in Figure 2-4 for particular values 
of p, mentioned above. The angle between the lines of marginal stability (corresponding 
to two roots z a , zj,) and the line Ai = is clearly determined by the fraction p = A1/A2. 
So by dividing the above equation with A2 and taking the limit Ai, A2 large we obtain: 

= -^P 5 (-! + ^) 3 ( 3 + V) + ^Pi- 1 + ^) 3 ( 3 + »)> ( 5 - 6 ) 

which has the real solutions p = ±1. So far away, the lines meet at an angle of tt/2. 
The same is the case in the k = 3 theory, where the curves of marginal stability (for the 
"basic" separatrix curve discussed in section 5.1) meet at an angle tt/2 at infinity. Now 
consider the Z 3 symmetric case as in Figure 2. For any pair of roots (z a , Zb) we have a basic 
separatrix curve joining them. Far away from the origin these curves meet at an angle of 
7r/2. Now, because of the Z 3 -symmetry, the angle between two neighboring curves must 
then be (7r/2)/3 = ir/6. Asymptotically we therefore have 12 domains. 

We start by counting the number of possible BPS states for the Z 3 -symmetric config- 
uration of roots, see Figure 2. Generally we will call z\ as root 1, Z2 as root 2 and so on. 
We start with the most symmetrical configuration, where the fourth root A is in the center 
of the triangle defined by the roots 1, 2 and 3. We call this small domain /. lis defined as 
the intersection of three domains: one where 1 is connected to 4 and 4 is connected to 3, 
but 1 and 3 is not connected (this comes from the basic separatrix curve connecting 1 and 
3), one where 2 is connected to 4 which is connected to 3, but 2 and 3 is not connected 
and finally one where 1 is connected to 4 and 4 is connected to 2 but 1 and 2 are not 
connected. This shows that the connectivity of the diagram in domain / must be of type 
B. The number of possible BPS states is therefore 3. The number of BPS states in the 
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other domains can now be determined by using the rules described in the last section in 
crossing the different separatrix curves. 

I— >II: In going to domain II one crosses the line -F143 = and since there was no 
connecting between 1 and 3 to start with these two roots gets connected by a BPS solution. 
The number of BPS states in 77 is then 4 and the connectivity is of type D. II— >-III: In 
going to domain III one crosses the line -F243 = and since there was no connecting 
between 2 and 3 to start with these two roots gets connected by a BPS solution. The 
number of BPS states in III is then 5 and the connectivity is of type E. Ill— >IV: In going 
to domain IV one crosses the line -F142 = and since there was no connecting between 1 
and 2 to start with 1 and 2 to will be connected such that all roots are connected and the 
number of BPS states is 6. The connectivity is of type F. IV— >V: In going to domain V 
one crosses the line .F314 = and since there was a connecting between 3 and 4 to start 
with, this domain wall disappears and instead 1 and 2 is connected. The connectivity is 
then of type E with 5 possible domain walls. V— >VI: In going to domain VI one crosses 
the line F i2 4 = and since there was a connecting between 1 and 4 to start with, this 
domain wall disappears. The connectivity is then of type D with 4 possible BPS states. By 
Z3-symmetry this analysis determines the possible number of BPS states in all domains 
and hence we have a complete determination of the possible domain walls and junctions 
for a potential with this particular symmetry. For this class of superpotentials, the number 
of BPS states varies from three to six. 

A similar analysis can be carried out for the Z2-symmetric case in Figure 3. When 
we simply plot the corresponding Eq. (|5.5| ) for all pairs of roots then we get more lines 
than is shown in Figure 3. However, some of these lines are lines of marginal stability, just 
like the real axis is for a quartic superpotential as discussed in section 5.1, and should be 
ignored. 
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Fig. 3: Z2-symmetric case. Connectivity of roots depending on the value of z± = A. The 
three fixed roots are located at z\ = —1,Z2 = 1 and Z3 = 3. 



However, here the three fixed roots are all colinear so the resulting diagram is very simple. 
To determine the possible BPS states in the different domains, one can start by taking 



1 < Z4 < 3 and real. In this case the configuration is known [ 3|J] : all roots are successively 
connected as shown in Figure 3. The configuration in other domains is then simply deter- 
mined by crossing the different curves of marginal stability. For this case the number of 
BPS states varies from three to five. 

The case including the Z4-symmetric potential is presented in Figure 4. At first glance 
this figure looks very complicated. However, it has some features common with Figure 2. 
For example, eight separatrix lines emanate from each critical point. For this choice of 
parameters, the number of BPS states varies from three (around the 'center') to six (at 
the Z4 symmetric point for example). So all possible connectivities are realized, except 
the minimal case of three BPS states (Figure 1-A) and Figure 1-C of course. 

We have found all the possible BPS states. Now let us focus on the junction configu- 
rations. As mentioned before, a triangle leads to a junction of three domain walls. If only 
one edge of the triangle is occupied, it is a 1/2 BPS state of a single domain wall. When all 
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three edges are occupied, then all the tensions will be balanced and this will lead to a 1/4 
BPS configuration of junctions of domain walls. More generally, one could have junctions 
of any number of domain walls. Here is how we can define a junction configuration in this 
case. First find the locations of the critical points in the IT-space. The integral curves will 
be straight lines between two critical points i,j in this space, and will have corresponding 
'soliton' number /i^, which can also be zero. Next pick any number of critical points in the 
IT-space, such that the successive connection of these points form a convex polygon. If all 
the edges have nonvanishing fj,ij, that is, if the polygon is closed then we have a nontrivial 
junction, and the inside of the polygon will have 1/4 supersymmetry. Each of the edges 
of the polygon will have 1/2 of the original supersymmetry, and only at the vertices, that 
is at the critical point, is all of the original supersymmetry preserved. Again, we see in 
this 'graphical' understanding of various SUSY configurations that there is no room for 



Fig. 4: Connectivity of roots depending on the value of Z4 = A. The three fixed roots are located 
at 01 = — 1, Z2 = 1 and Z3 = — 1 + 2i. In the empty domains the connectivity is of type E. 

If any of the fiij along the edge of the polygon is zero, then we cannot define the 'inside' 
of the polygon and there will be no junction. We will have just domain walls with extend 



3/4 BPS states in the W-Z model 
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to infinity (in the coordinate space) and which never join. The number of preserved 
supercharges will be two. 



6. Discussion 

So far we have been discussing the possible BPS states and junctions in the W-Z model. 
We have summarized our result in Figures 2-4 where we can easily read off the number 
of BPS states as well as possible BPS junction configurations for a given deformation 
parameter. So what is the use of all this? First of all, we have used a method general 
enough to be utilized for counting BPS states for other types of superpotentials. Secondly, 
the BPS data of W-Z models (or those with other superpotentials) which can be obtained 
from higher dimensional theories will reflect the BPS data of the original theories. 

Apart from these practical things, we would also like to point out some of the pos- 
sible connections to works done in the context of string compactifications and also brane 
configurations. Due to the relation to Calabi-Yau compactifications we can reinterpret our 
results as that of counting numbers of BPS D-branes wrapped around supersymmetric 
cycles. On top of each domain wall there is a 'sphere' wrapping around a supersymmetric 
cycle, whose radius vanishes at the critical points. This is reminiscent of toric geometry: 
We have vanishing spheres at the critical points and have finite radius cycles over the line 
interpolating two critical points. That is spheres in the internal dimension over the domain 
walls will be revealing some of the structures of Calabi-Yau spaces. In particular, it has 
been shown that certain toric geometries, which has vanishing cycles, can be translated 



into a brane configuration |35| . Thus another very interesting application comes from the 
T-dual picture of the Calabi-Yau compactifications, i.e. the brane configurations. As an 
example consider the following situation. The brane configuration for M = 1 SU(N C ) su- 
persymmetric YM is given in Type IIA string theory with N c D4 branes extending between 
two sets of coincident NS5 branes as follows. With the coordinates 

s = x 6 + ix 10 , v = x 4 + ix 5 , w = x 8 + ix 9 . (6-1) 

the D4 brane is located at v = w = and extended in the s-direction. The NS5 brane 
is located at s = w = and is extended in the ^-direction, and the NS5'-brane is at 
v = 0, s = L and is extended along the w-direction. Consider a configuration of m 
coincident NS5 branes connected by N c D4 branes to m' coincident NS5' branes. There 
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will be two adjoint superfields $, which describe fluctuations of the fourbranes in the 
w and v directions respectively, whose classical superpotential is 



a „ L1 a' 



W = — Tr$ m+1 + Tr$ m ' +1 + • • • . (6.2) 

m + 1 m' + 1 

Imagine having the k NS5 branes in the (x 8 ,x 9 ) plane at k different points Wj. Since 
the {wj} correspond to locations of heavy objects they appear as parameters rather than 
moduli in the gauge theory description and give rise to a polynomial superpotential for $ 
where W'(3>) = aYYjLii^ ~ w i)- This shows how superpotentials of the form discussed in 
this paper can arise from brane configurations. 

Another very interesting result can be obtained with similar methods in the study of 



BPS states of Argyres and Douglas superconformal theories Pq] ||37|| [p8| , as in Ref. |34 
In fact, if we consider a degenerate choice of polynomial, where P m = (dW/dx) 2 , the 
problem becomes identical to the problems we have discussed here. Exact equations for 
the separatrix curves can be obtained but will be quite complicated and involve certain 
elliptic functions. 

As discussed in section 3, when we consider Type IIA string theory compactified on a 
Calabi-Yau fourfold we obtain a 1+1 dimensional effective theory which gives the vacuum 
structure and the D4 branes wrapping around the supersymmetric cycles give solitons 
interpolating the vacua. If we start with M-theory, which is the strong coupling regime 
of the Type IIA theory, we end up with an effective 2+1 dimensional theory, with similar 
vacuum and domain wall structure. However, there is something more. Due to one more 
space dimension, the vacua can arrange such that there can be junctions of the domain 
walls. From the point of view of string theory this extra dimension is a nonperturbative 
effect. Thus having a full understanding of lower-dimensional integrable models might not 
guarantee an understanding of higher-dimensional integrable model, just as understanding 
fully perturbative field theory does not guarantee any insight into a fully nonperturbative 
field theory. 

The superpotential we have studied in this paper is the simplest kind involving only 
one type of field. There are many extensions that can be made with multiple species of 
fields. One nice extension would be the study of the D — E series [[58] of singularities and 
the corresponding W-Z models. In the case of W-Z models of A n type, one always has a 
single type of domain walls between two critical points, because there is only one type of 
complex scalar field in the theory. However, if we have multiple species of scalar fields we 
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have the possibility of multiple types of domain walls between the critical points. It would 
be interesting to generalize the method used here to study these systems and also find 
junctions of multiple species of branes. Theories such as the CP n models have multiple 
species of domain walls between critical points, which can be labeled by a group theory 
index. So when we consider junctions of a multiple of these domain walls, perhaps only a 
certain combinations will lead to a BPS junction. This certainly deserves a further study. 

There are still some open questions, we would like to answer in the near future: How 
do we describe junctions of domain walls in the higher-dimensional Calabi-Yau geometry? 
Are stable junctions classified by some topological class, related to the higher-dimensional 
geometry? 
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